Coulomb gap in one-dimensional disordered electron systems 
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The density of states of one-dimensional disordered electron systems with long range Coulomb 
interaction is studied in the weak pinning limit. The density of states is found to follow a power 
law with an exponent determined by localization length, and this power law behavior is consistent 
with the existing numerical results. 
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/jTiroducfion-Recently there have been much interests 
in the one-dimensional (ID) electron systems motivated 
by the development of carbon nanotube technology.!!! 
In ID electron systems electron-electron interactions 
play very important roles, leading ta phases different 
from the conventional Fermi liquids .□ Repulsive short 
range electron-electron interactions cause Luttinger liq- 
uids (LL)J3 while long range Coulomb-interaction (LRCI) 
is believed to cause a Wigner crystal. B The Wigner crys- 
tal phase of spinful electrons is characterized by quasi 
long range order of 4kp charge density components. 
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Eq. ([!]) should be compared with the power law depen- 
dence x~ a of LL. Eq. (0) also indicates the suppressed 
quantum fluctuations of charge densities, and this feature 
manifest in the density of states (DOS) D(uj)\ 
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where to is measured from Fermi energy, and E c is a cut- 
off energy. Note that DOS of Eq. (||) decays faster than 
any other power law. For LL, DOS follows a power law 
|w| 7 with a nonuniversal positive 7. The physics of ID 
Wigner crystaLis analogous to that of charge density 
wave (CDW)eI'EIeI apart from the presence of quantum 
fluctuations and LRCI. 

Impurities, either a few or many, change the physi- 
cal properties of ID electron systems qualitatively. For 
a single impurity in an electron system with repulsive 
interaction, the back scattering of electrons with the im- 
purity becomes strong at low energy, and it effectively 
divides the system into two pieces. I3fcl 

For noninteracting electrons in random disorder (many 
impurities), all the states are known to be localized due 
to repeated back scattering.^ DOS for a disorder with a 
Gaussian .distribution (with zero mean) can be calculated 
exactlyllEEd 
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In the strongly localized electron systems with LRCI, 
where the overlap of wave functions can be neglected, 
electrons can be treated classically and DOS exhibits a 
Coulomb gap of the form:E3 



(4) 



We note that Eq. (Q) has been derived under the assump- 
tion that the localization length is much smaller than 
the interparticle distance. Both disorder and LRCI push 
a ID electron system to classical limit but in a different 
manner as is reflected in the form of DOS Eq. (g) and Eq. 
(||). In this paper we report a result on the DOS at low 
energy for the ID disordered electron system with LRCI 
when the localization length is larger than the interparti- 
cle distance or interimpurity distance, [more precisely the 
weak pinning limit, see below ] Following the analys es . . on 
a pinned Wigner crystal by Maurey and GiamarchiESEJ, 
employing a simplified model, and using semiclassical ap- 
proximation we find that DOS follows a power law at low 
energy. 
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where the exponent r/ is basically determined by the lo- 
calization length, [see Eq. Jig)] This power law behavior 
is consistent with the existing numerical results JIj Eq. 
(||) is the main result of this paper. 

Model- We consider a spinless electron system for sim- 
plicity. Such a system can be realized in organic chains 
and quasi-lD quantum wires in strong magnetic field.t£l 
The Hamiltonian consists of three parts. 
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The operator iPr(iPl) is the right-moving (left-moving) 
electron operator. The continuum chiral electrons and 
lattice electron operators are related by 
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where a is lattice constant. 

Pr —'■ yr^r '■ is the normal ordered right-moving 
edge electron density operator (pl is similarly defined), 



and p(x) = p R (x) + p L (x). V(x) = %- v J +d2 is the 
Coulomb interaction, d is the transverse size of quantum 
wire which we take to be the same a for simplicity, e 
is a dielectric constant. The Coulomb matrix element is 
V{k) = *fK (a\k\) ~ ^ln^ for \k\a > 1 and K is 
the modified Bessel function. The impurity Hamiltonian 
is given by 



^imp-E^)^)^ 



where 7 is a numerical constant characterizing short- 
range interaction and quantum fluctuations and is 
impurity concentration. The logarithmic factor is due 
to LRCI, which enhances the pinning length. The en- 
hanced pinning length implies the more rigid system and 
the more difficult pinning by impurties. The system be- 
comes much more ordered and the fluctuations around 
the ground state is much less important. The expression 
for pinning length Eq. (13) has been derived under the 
assumption of weak pinning Lq Ti^ . 

Beyond the pinning length scale the phase coherence of 
9 is lost, thus the correlation function of CDW operator 
is expected to decay exponentiallyQ (\x — x'\ > Lq) 



dxWi{x)[p(x) + e 2lkFX i)l(x)ip R {x) +H.c],(9) (cos[20(x, t)] cos[20(a/, t)]) ~ e -\*-*\/L°, (14) 



where Wi{x) is the impurity potential. The first term 
in the bracket of Eq. (g) is the forward scattering term 
which can be neglected compared to back scattering at 
low energy. The impurity potential Wi(x) is chosen to 
be 
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where Xj's are the random locations of impurities. The 
interactingpdectron systems can be bosonized in a stan- 
dard way.aliU The phase fields and bosonization formulas 
are given by 
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The bosonized Hamiltonian of the system is 
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where po is the average density of electrons. El 

Pinning Length-The pinning length Lq is a length scale 
over which the phase field 9{x) in the ground state varies 
with 59 7r ip order to take advantage of the impurity 
potential. EJ'til The pinning length corresponds to the lo- 
calization length of the electron system. The pinning 
length can be obtained by maximizing the energy, gains 
from impurity, elastic, and Coulomb energies .EMa In- 
cluding the effect of quantum fluctuations using the self- 
consistent harmonic approximation, the pinning length 
is given by 
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Eq. ([IT]) implies that the system breaks into segments 
whose typical length is given by the pinning length Lq. 
Thus at low energy we can consider a typical segment 
with length Lq and calculate DOS averaged over the seg- 
ment. Since the tunneling between segments is strongly 
suppressed at low energy we have to fix the value of 9 at 
the boundary of a particular segment J3(Dirichlet bound- 
ary condition) 

Approximate Low Energy Model- Let us construct a 
model for a segment which is valid at low energy. First of 
all the logarithmic divergence of Coulomb matrix element 
V(k) is cut by k ~ l/L . Then the Coulomb energy term 
can be expressed as 
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The effective Hamiltonian of a segment becomes 
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where V(9, < Xj < L ) is some potential energy term 
which is optimized by a certain random average phase 
value 9 pt and it depends on 9 only through e ±2l9 .[see 



the last line of Eq. fll2[)l Note that 9 — * 9 + ir is then a 
symmetry of the system. In the weak pinning limit we 
are considering the phase 9 varies rather smoothly in the 
range \9 — 9 opt \ < ir at low energy over a length scale Lq. 

Calculation of DOS- The electron Green function in 
imaginary time is defined by 

GR/L{x,y,n - t 2 ,Wi) = -(T r i/j R/L (x,T 1 )ip R/L (y,T 2 )) 

for a particular realization of impurities Wi . After aver- 
aging ove impurities 

GR/ L (x-y,T 1 -T 2 ) = (G R / L (x,y,T 1 -T 2 ,W I )) Wl - (18) 
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DOS is defined by 



D(u) = —Tm / dre l ^G R/L (0,T) 
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The above expressions are formally exact but difficult 
to calculate. Here instead we will calculate the electron 
Green function at a point y in the segment (away from 
the boundary) using the Hamiltonian Eq. (16). 



(Gn/L(y,y,Ti - t 2 )) Wi - 
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For the computation of electron Green function, the La- 
grangian formulation is more convenient. 



in 







S scs = I dx j dr-d T 9d x 4>+ \ drH scg . (21) 



Using the bosonization formula Gn/L^y, t% — T2) can 
be expressed as 



G R (y,Ti - r 2 ) = 



dxdr 



Sex 



J D[6,(j)}e- S ^- 
9(x,t) + (/>(x,t) J(x,t), 



J(x, r) = S(x - y) 6{t - n) - 5(t - r 2 ) . (22) 



In Eq. (|22|) the dual phase field <j) can be integrated out 
explicitly. [Remember that the potential term V(9,Xj) 
depends only on 9 not on </>.] The electron Green function 
becomes 



Gn(y, n - t 2 ) 
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Because our system is in the deep classical regime the 
dominant low energy processes would come from the 
quantum tunneling between classical vacua. Remem- 
bering that the phase field 9 is an angular variable, the 



classical vacua are charerized by # vac ~ {6> opt + nn, n — 
integer}. This is because the potential V(9,Xj) depends 
on 9 only through the form e ±2 * e , so that it possesses the 
symmetry 

9^9 + n (27) 

Since the phase field 9 varies on the order ir in each 
segment one can expect that the most domiant pro- 
cess would be the quantum tunneling between ^opt and 
9 op t ± 7r vacua. The quantum tunneling processes be- 
tween classical vacua can be described by the solutions 
of classical equations of motion in imaginary time, which 
are also called instantons\}3 

Now we have to find the classical (imaginary) time 
dependent solution which minimizes So + Sje . Since Sjg 
is linear in 9 it plays a role of the external source field. 
We argue below that the first term of Sjg [ Eq. (|25|) ] is 
the source for the vortex and anti- vortex configuration of 
the (angular) phase field 9. For the moment let us neglect 
the potential term V(9, Xj). Then the phase field 9 and 
its dual phase field <fi are related by the Cauchy-Riemann 
equatiorilj [we set u — 1 for simplicity from now on] 



d^<j) — -iefj.udvB, fi,v = T,x, 



(28) 



where e M „ is a totally antisymmetric tensor. The first 
term Eq. ( p5| ) stems from the source of dual field 4>.[see 
Eq. (H|)] Including only one source at (ri,y) the cj> field 
satisfies the classical equation of motion. 

-(d* + d%)(i>(x,T) = in6(x-y)8(T-T l ). (29) 

Combining Eq. ( f28|) and Eq. (|2|) we deduce 

nS(x - y)S{r - n) = e liy d tl d l ,9{x, t), (30) 

Naively the right hand side of Eq. ([30]) vanishes but 
this is not necessarily true for the topologically singular 
configuration such as vortex. The vortex nature of 9 is 
demonstrated by the calculation of vorticity using Eq. 
(p0|) and Stokes theorem. 



dx^ 8^6 



d 2 xeav dnd„9 = l-2ir 



(31) 



which characterizes the vortex nature of (9 .£3 The unusual 
factor of 1/2 is due to the aforementioned symmetry 
9 — * 9 + 7r. Remembering that we have another souce at 
(T2, y) with a charge opposite to that at (ri, y)[see the last 
line of Eq. (p^)] , we conclude that the vortex-antivortex 
configuration of 9 phase field will dominate the electron 
Green function at long time in the classical regime. Be- 
fore proceeding to the explicit vortex-antivortex instan- 
ton solution, let us discuss another classical field configu- 
rations generated by the second source term of Eq. (|2q). 
It satisfies the standard Laplace equation in two dimen- 
sion with delta function sources, and its explicit form is 
given by 



%( T ,x)~-ln (x 



K {x _ y) 2 + {T _ Tl f 
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In the configuration of Eq. (32) the range of variation of 
phase field 9 is much large than 7r for long time limit \t\ — 
t 2 | — > oo. Thus the configuration 8 (x) is suppressed by 
the potential term V(6,Xj), and it also does not exhibit 
the angular nature of phase fields. This is an analogue of 
the spin- wave degrees of freedom of XY model Ji3 Thus, 
the configuration of Eq. ( |3^ ) does not contribute the 
electron Green function in long time limit. 

The explicit vortex-antivortex solution in complex co- 
ordinate is given byliij 



ai 29( Z ) _ z - (n + iy) 



Z = T + IX 



z- (t 2 + iy) 

Substituting Eq. (|33|) into the action Sg we obtain 



(33) 
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where To is the short time cut-off provided by the con- 
tributions from the vortex cores. Eq. (Em|) implies Eq. 

®- _ 

To complete our analysis let us compute Sj Eq (26). 
The divergent k summation at small k should be cut 
by I/Lq. The frequency summation diverges at large 
frequency (short tigie) which is cut by some high energy 
scale such as vpn^a Thus, Sj does not contribute to the 
electron Green function in long time limit. 

According to Eq. (j^) the exponent of DOS is larger for 
the longer localization length, which-jjs consistent with 
the numerical results by Jeon et al.£3 The exponent r\ 
can be expressed as 



1 2c 1 , Lo 
V = T^z x — x - x In— , 
137 vf e a 



(35) 



where c is the speed of light. In the weak pinning limit 
it is reasonable to take In — to be 4-6. Choosing a typ- 



ical Fermi velocity near 10 7 cm/sec and e ~ 1 — 5 for 
a quantum wire, we can estimate the exponent of DOS 
yT + rj/2 to be around 3-6 which is also consistent_ssiith 
the exponents obtained from the numerical studies .113 

More precisely the Eq. should be averaged over 

the impurity configurations in the segment. Inside the 
segment the phase field varies smoothly in space and the 
impurity average is not expected to bring in any singular 
effects which would invalidate Eq. (Q) , since no infrared 
divergence can occur through impurity averages for a fi- 
nite segment. 

Sufficiently strong short range interaction can also 
push the system to classical regime (K << I), where the 
result Eq. (0) is applicable. On the contrary, the result 
Eq. (H) is not applicable to the case of noninteracting 
disordered electrons (K = I) since the condition of the 
classical limit (K << 1) is not met. Indeed the DOS of 
noninteracting disordered electrons is finite at the Fermi 
level (See Eq. (§)) contrary to Eq. (|). 

Summary- We have studied analytically the DOS of 
the disordered ID electron system interacting via long 
range Coulomb interaction employing a simplified model 
and semiclassical approximation. The DOS is found to 
follow a power law with a nonuniversal exponent which is 
basically determined by localization length. The power 
law dependence is also consistent with the existing nu- 
merical result. 
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